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ABSTRACT 

We perform test-particle simulations in a 2D, differentially rotating stellar disk, subjected to a 
two-armed steady state spiral density wave perturbation in order to estimate the influence of spiral 
structure on the local velocity field. By using Levenberg-Marquardt least-squares fit we decompose the 
local velocity field (as a result of our simulations) into Fourier components to fourth order. Thus we 
obtain simulated measurements of the Oort constants, A, B, C, and K. We get relations between the 
Fourier coefficients and some galactic parameters, such as the phase angle of the Solar neighborhood 
and the spiral pattern speed. We show that systematic errors due to the presence of spiral structure 
are likely to affect the measurements of the Oort constants. Moderate strength spiral structure causes 
errors of order 5 km/s/kpc in A and B. We find variations of the Fourier coefficients with velocity 
dispersion, pattern speed, and sample depth. For a sample at an average heliocentric distance of 0.8 
kpc we can summarize our findings as follows: (i) if our location in the Galaxy is near corotation 
then we expect a vanishing value for C for all phase angles; (ii) for a hot disk, spiral structure 
induced errors for all Oort constants vanish at, and just inward of the corotation radius; (iii) as one 
approaches the 4:1 Lindblad resonances |C| increases and so does its variation with galactic azimuth; 
(iv) for all simulations |C|, on average, is larger for lower stellar velocity dispersions, contrary to recent 
measurements . 
Subject headings: 



1. INTRODUCTION 

Jan Oort proposed (^Oortlll927a^ and measured (jOortI 
ll927bD the now called Oort constants A and B, over 
seven decades ago by assuming that stars in the Galaxy 
moved on circular orbits. Later on his ana lysis was gen- 
eraliz ed to the non-axisymmetric case (e.g.. lOgorodmkovl 
119321 ) resulting in the definition of two more constants: 
C and K. The importance of the Oort constants (OC) is 
in their simple relation (at least in the case of vanishing 
random motions) to the galactic potential. 

Many investigations have been conducted in the at- 
tempt to measure the OC. Howeve r, not all the results 
are in agreement (for a review see iKerr fc Lvnden-Belll 
1198 6^. which is mainly due to the absence of a complete 
proper motions catalog. Recent studies have failed to 
improve measurements of the OC. 

Another reason for the diversity am ong the values de- 
rived for the OC has been discussed bv lOUing &: DehnenI 
(|2003f l. The authors present an effect which arises from 
the longitudinal variations of the mean stellar parallax 
caused by intrinsic density inhomogeneities. Together 
with the reflex of the solar motion these variations create 
contributions to the longitudinal proper motions which 
are indistinguishable from the OC at < 20% of their am- 
phtude. OUing & Dehncn (2003) (hereafter O&D) mea- 
sured the OC using proper motions from ACT/Tycho- 
2 catalog and corrected for the "mode mixing" effect 
described above, using the latitudinal proper motions. 
They found that the OC vary approximately linearly 
with asymmetric drift (or age) from the young-stars val- 
ues to those appropriate for old populations. This trend 
may result from perturbations induced by the Galactic 
bar or spiral structure. In this paper we investigate the 
latter possibility. 



Deviations from the axisymmetric values of the OC 
have not received a satisfactory interpretation although 
they are often assumed to arise from a non-axisymmetric 
Galactic potential perturbation such as spiral structure 
or/and the Galactic bar. 

Given their preeminent importance, we feel that a bet- 
ter understanding of the effect of spiral perturbation on 
the measurements of the OC is necessary. Spiral arms 
perturb the local velocity field corrupting local measure- 
ments of the OC. Thus we look for the effect of spiral 
structure on the OC A, B, C, and K and their deviation 
from the axisymmetric limit. We hope that this would 
give us some constraints on the spiral structure. Our ap- 
proach applies to the kinematics of the Local Standard of 
Rest (LSR) , while at the same time it constitutes a gen- 
eral analysis of a spiral density wave perturbing a stellar 
disk. 

Our simulation model is described in ij2l We present 
a derivation of the OC in iJ3]and describe our numerical 
measurement technique in f|4l Results are presented in 
§35HZl Discussion of observational measurements of the 
OC is given in ^ and summary follows in ^ 

2. NUMERICAL MODEL 

2.1. Notation and units used 

For simplicity we work in units in which the galactocen- 
tric distance to the annulus in which stars are distributed 
is To = 1, the circular velocity at tq is = Ij and the 
angular velocity of stars is flo = vq/tq — 1. Because we 
assume a flat rotation curve throughout this paper the 
initial circular velocity is = 1 everywhere. One orbital 
period is 27rro/rio — 27r. In our code time is in units of 
1/0. The velocity vector of a star is (u, w), where m, v are 
the radial and tangential velocities in a reference frame 
rotating with vg. Consequently, the tangential velocity 
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of a star in an inertial reference frame is = vq + v. We 
refer to the radial and tangential velocity dispersions as 
au and ay, which are defined as the standard deviations 
of u and v, respectively. 

The heliocentric distance, azimuth, radial and trans- 
verse velocities are denoted by d,l,Vd, and vi, respec- 
tively. As it is commonly accepted, we also refer to / as 
the Galactic longitude. We write the longitudinal proper 
motion as /i;. 

The azimuthal wavenumber of the spiral wave, m, is an 
integer corresponding to the number of arms. The grav- 
itational potential perturbation amplitude of the spiral 
density wave is denoted as and the pattern speed is 
fig- The parameter a is related to the pitch angle of the 
spirals, p, with a = mcot(p). a is negative for trailing 
spirals with rotation counterclockwise. In this paper we 
only consider trailing spiral density waves. 

Throughout this paper we refer to the dynamical Local 
Standard of Rest as LSR, which is a point at the position 
of the Sun moving with a circular velocity vq. The Sun's 
velocity vector used in the derivation of the OC is 
defined as (u,v^) = {Uq,Vq). 

2.2. Equations of motion 

We consider the 2D motion of a test particle in the 
mid-plane of a galaxy in an inertial reference frame. In 
plane polar coordinates (r, (p) the Hamiltonian of a star 
can be written as 

H{r,(j),pr,P4>,t) = Ho{r,pr,p^) + $i(r,0,i) (1) 

The first term on the right hand side is the unperturbed 
axisymmetric Hamiltonian 

Hoir,Pr,p^)= ^ + ^ + ci>o(r) (2) 

where we assume the axisymmetric background potential 
due to the disk and halo has the form ^o{t) = Vq log(r), 
corresponding to a flat rotation curve. 

The imposed gravitational potential perturbation, 
$i(r, (j), t) is due to a spiral density wave and is discussed 
in the next section. 

The equations of motion derived from the Hamiltonian 
(eq. [1]) are integrated particle by particle. 

2.3. Perturbation from spiral structure 

In accordance with the Lin-Shu hypothesis (|Lin et al.l 
Il969t l we treat the spiral pattern as a small logarithmic 
perturbation to the axisymmetric model of the galaxy 
by viewing it as a quasi-steady density wave. We expand 
the spiral wave gravitational potential perturbation in 
Fourier components 



$i(r,0,i) 



^i[a In r—m{(l) — 0,st — (l)m)] 



(3) 



and its corresponding surface density 

We assume that the amplitudes, £,„,!],„, and the pitch 
angle are nearly constant with radius. The strongest 
term for a two-armed spiral is the m = 2 term. Thus 
only the terms corresponding to m — 2 are retained. 



Upon taking the real part of eq. [3] the perturbation due 
to a spiral density wave becomes 

<I>i(r, <j), t) — €s cos[a Inr — to((/) — VLgt — (t>o)\- (5) 

The direction of rotation is with increasing </> and a < 
ensure that each patt ern is trailing . 

The recent study of lValled (|2005f l provides a good sum- 
mary of the many studies which have used observations 
to map the Milky Way disk. Ccpheid, HI, CO and far- 
infrared observations suggest that the Milky Way disk 
contains a four-arm e d tigh tly wound structure, whereas 
iDrimmel &: Spergell (|2001h have shown that the near- 
infrared observations are consistent with a dominant two- 
armed structure. A dominant two-armed and weaker 
four-armed structure, both moving at the same pat- 
tern sp eed, was previously proposed by Amaral & Lepini 
()1997f ). A combination of primary two- and weaker 
four-armed spiral structur es, moving at different patt ern 
speeds was considered by iMinchev fc Quillenl (|2006f ) as 
a possible heating mechanism. Here we consider a two- 
armed spiral density wave only and differ the four-armed 
case and the two+four arms combination to future inves- 
tigations. 

In general, an individual stellar orbit is affected by 
the mass distribution of the entire galaxy. However, if 
tight winding of spiral arms is assu med only l ocal g ravi- 
tational forces need be considered. iMa et ahl (|2000 '1 find 
kr for several spiral galaxies of various Hubble types to 
be > 6, thus the tight- winding, or WKB approximation, 
is often appropriate. In the above expression k is the 
wave- vector and r is the radial distance from the Galac- 
tic Center (GC), related to the pitch angle p through 
cot(p) = \kr/m\. This gives a in terms of fc as a ~ |fcr|. 
In the WKB approximation the amplitude of the po- 
tential perturbation Fourier component is related to the 
density perturbations in the following way 

-27rGS]sro 
« p-r^ (6) 

(|Binnev fc Tremaind Il987f ) . The above equation is in 



units of 



Es is the amplitude of the mass surface 



density of the two-armed spiral structure. 

2.4. Resonances 

Of particular interest to us are the values of fig which 
place the spiral waves near resonances. The Corotation 
Resonances (CR) occurs when the angular rotation rate 
of stars equals that of the spiral pattern. Lindblad Res- 
onances (LRs) occur when the frequency at which a star 
feels the force due to a spiral arm coincides with the 
star's epicyclic frequency, k. As one moves inward or 
outward from the corotation circle the relative frequency 
at which a star encounters a spiral arm increases. There 
are two values of r for which this frequency is the same 
as the epicyclic frequency and this is where the Outer 
Lindblad Resonance (OLR) and the Inner Lindblad Res- 
onance (ILR) are located. 

Quantitatively, LRs occur when f2s = f2o i n/m,. The 
negative sign corresponds to the ILR and the positive - 
to the OLR. Specifically, assuming a flat rotation curve, 
for the 4:1 ILR = 0.65f7o, for the 4:1 OLR = 
l.SSf^o, for the 2:1 ILR 0.^ = 0.317o, and for the 2:1 OLR 
ils = 1.7r2o. Corotation of each spiral pattern occurs at 
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2.5. Simulation method 

To investigate the effect of a spiral density wave on the 
local velocity field we need to simulate a large number of 
stars. In this paper we consider an initially cold or hot 
stellar disks as described in ^'2.6\ We distribute parti- 
cles in the annulus (0.6ro, 1.4ro), when simulating a cold 
disk, or (0.4ro, 1.4ro) for the hot-disk case. We reduce 
the computation time by recording the positions and ve- 
locities of test particles 33 times per orbit for three orbits 
after the initial seven rotation periods. Position and ve- 
locity vectors are saved only if particles happened to be 
in the statistically important to us region (O.Srg, 1.2ro). 
We require the final number of stars in the case of the 
cold disk to be 5 X 10^. For the hot-disk simulation more 
particles are needed to beat the random motions: we 
require the output to be 1.5 x 10^. After utilizing the 
m-fold symmetry of our galaxy, we end up with a work- 
ing sample of 5 x 10^ (cold disk) or 1.5 x 10^ (hot disk) 
stars in our fictitious Solar Neighborhood (SN). 

This way of creating additional particles is justified by 
the fact that a steady state (constant pitch angle and 
angular velocity) spiral density wave does not change 
the kinematics of star s once it has been gro wn. As was 
demonstrated by .Minchev fc Quillenl ()2006f ) test parti- 
cles remain on the approximately closed orbits assumed 
during the growth of the spiral pattern unless the pattern 
speed places stars at an exact Lindblad resonance. 

2.6. Initial conditions 

The density distribution is exponential, ^(r) ^ e"''/'"'', 
with a scale length rp = 0.3 8ro, c onsist ent with the 
2MASS photometry study bv lOihal ()200H ). We would 
like to see how spiral structure changes the velocity field 
in the case of an initially cold or initially hot stellar disk. 
To simulate the cold disk we send stars on circular orbits 
with the same initial azimuthal velocity Vtf, = vq + v = vq 
consistent with a flat rotation curve. For the hot disk we 
give an initial radial velocity dispersion in the form of a 
Gaussian distribution. The Milky Way disk is known to 
have a radial velocity dispersion which decreases roughly 



exponentially outwards: 



~ e 



-r/r^2 _ jjj accordance 



with this we implement an exponential decrease in the 
standard deviation of the radial velocity dispersion of 
stars with radius, with a scale length r„2 ~ 0.45ro 
(jLewis k Freeman! 1 1 9891 ) . We set (t„ = 40 km/s at tq. 

The amplitude of the spiral density wave, Cs, is ini- 
tially zero, grows linearly with time at < t < ti and 
is kept constant after i = <i = 4 rotation periods. Thus 
the spiral strength is a continuous function of t, insur- 
ing a smooth transition from the axisymmetric to the 
perturbed state. 

Table [T] summarizes the parameters for which simula- 
tions were run. 

2.7. Numerical accuracy 

Calculations were performed in double precision. We 
checked how energy was conserved in an unperturbed 
test run. The initial energy was compared to that cal- 
culated after 40 periods and the relative error was found 
to be \AE/E{0)\ < 8.48 x 10"". We also checked the 
conservation of the Jacobi integral J in the presence of 
spiral structure. The relative error found in this case was 
|AJ/J(0)| < 3.11 x 10-12. 



3. DERIVING THE OORT CONSTANTS 

Consider a star moving in the gravitational potential 
described in H2.2[ We assume that its velocity vector with 
respect to an inertial frame is {u,v^) where u is radial 
and Vci, = vq + V is tangential. Let u be positive toward 
the GC. We define to be positive in the same direction 
as the galactic rotation. From the geometry of Fig. [T]we 
can write the radial and transverse velocities of this star, 
Vd and vi respectively, as seen by an observer sitting on 
the Sun, in terms of its galactocentric velocities (m, v^) 
as 

vi{l) ~ Uq sin I — Vq cos / -I- cos(Z + 0) — m sin(/ + 0) 
VdiJ-) = —Uq cos / -I- Vq sin / + sin(/ + (!))+ ucos(/ -I- (j)!)) 

where / is the heliocentric azimuth, r is the radius of the 
star from the Galactic Center (GC), and (p is the angle 
between the star/GC and the Sun/GC vectors. {Uq, Vq) 
is the velocity vector of the Sun. Angle definitions lead 
to 



sin(/ + </)) =r sint/) 
d cos / -|- r cos 4> — ro. 



(8) 



Using the laws of sines and cosines together with eq. ([5]) 
we eliminate (j) from eq. ([7]) in favor of d, the heliocentric 
distance of the star, arriving at 

vi{l) — Uq sin / — Vq cos / -I — -{ro cos / — d) ro sin/ 

r r 

Vd(l) — "UQsml — VqCOsI H — (rocos/ — d) H rosin/ 

r r 

To first order in d or (h. 



and infinitesimals 



dr : 



dcosl 



-dcosl 



d(/)— — sin/. 
ro 

We now Taylor-expand and u to first order about the 
Sun, (ro,0o): 

V4,{r, (/)) =-y^(ro,0o) + ^(roiMdr + ^(ro,0o)rf0 

du du 
u{r, 4>)=u{rQ,(j)Q) -f —{rQ,(l)Q)dr + —{rQ,(j)Q)d(j). 

Assuming circular motion, the motion of the Sun implies 
that Vff,{ro,(j)o) = Vq and u(ro,(/)o) — Uq. Using the 
above expressions for v, r, dr, dip to first order in d we 
can write the velocities vi , as 

vi{l) = Uq sin(/) - Vq cos(/) -I- ^dcos(2/) - Cdsin(2/) -I- Bd 
vd{l) = -Vq sin(/) - Uq cos(/) + Adsm{2l) + Cdcos{2l) + ^ 
where 

1 du 

r dr r dtp 
Vfjj dv^ 1 du 
r dr r dcp 
u du 1 dVrf, 
r dr r dcp 
u du I dv^ 
r dr r dcp 



2A^ + ^-± 



2B\ 



2Ce 



2K\ 



(10) 
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These equations are consistent with eq. (3) of O&D in 
slightly different notation. 

We see from eqs. IHlthat A and C can be determined 
from either the longitudinal proper motions /z/ = vi/d or 
radial velocities, Vd- When A and C are estimated from 
radial velocities their values are sensitive to the errors in 
the distances, whereas proper motions provide a distance 
independent way of measuring them. On the other hand, 
B can only be measured from proper motions and K 
only from radial velocities. It is interesting to note that 
due to the uncert ainty in the definition of a n inertial 
coordinate system (jCuddeford fc Binnevlll99J . O&D) B 
is often estimated ^ from A and the Oort ratio a^"^ ja^ 

^ = -^ (11) 
a„2 A-B- ^''^ 

If there is no streaming then = §7 = = 0, so 
that C = = 0, and A, B reduce to the well known 
axisymmetric Oort constants 

dr 



2A, 



ro 



2B, 



dr 



(12) 



If, in addition, the axisymmetric stellar disk is cold and 
the rotation curve is flat we can define 



ro 
ro ' 



2Bo. 



(13) 



As seen in the derivation above, the OC are not really 
constant unless they are measured in the SN. In fact, 
they may vary with the position in the Galaxy (r, (p) . The 
reason for variations of the OC with Galactic radius, if we 
assumed the Galaxy were axisymmetric, is the fact that 
the number density and velocity dispersions are functions 
of r. In the case of a non-axisymmetric Galaxy, due to 
spiral arms, bar, etc., we would expect a variation with (j) 
as well. Thus, the Oort constants have often been called 
the Oort functions. 

4. NUMERICAL MEASUREMENT OF THE OORT 
CONSTANTS 

First we move to a reference frame centered on (ro, 0o) 
where is the phase angle of the SN defined in fig. 
[21 On the circle of radius rp we change the position of 
the SN between the convex and the concave spirals, thus 
treating 0o as a free parameter. We then select a sam- 
ple of stars at a particular distance and annulus width 
from the point (rp, (t>o) which we identify with the LSR. 
The radial and transverse velocities of the SN stars, Vd 
and vi respectively, as seen from the LSR are written in 
terms of the galactocentric velocities {u,Vif,) (cf. eq. [9]). 
Each annulus of stars is divided into 30 "heliocentric" 
azimuthal bins with the Galactic longitude I — point- 
ing toward the Galactic center. We obtain the functions 
fLi = vi{l)/d{l) and Vd{l)/d{l) where /x; is the longitudi- 
nal proper motion and the bars indicate average values. 
Hereafter we omit the bars. 

^ The relation between the Oort ratio and the Oort constants 
is only justified in the case of well mixed, low velocity dispersion 
populations. Third order m oments may contribute sign ificantly if 
old stellar samples are used llCuddeford fc BinneyjUoM) . 



4.1. Fourier expansion of fii andvd/d 
We decompose /j,; = vi/d and Vd/d into Fourier terms: 

n 

fii = coi + ^(sii sin?; + cu cos il). 
1=1 

n 

— ^cod + y^{sidsmid + cidcosid). (14) 
d ^ 

By using Levenberg-Marquardt least-squares fit we ob- 
tain the expansion coefficients up to n = 4. The OC A 
is given by C21 and S2d'i C is given by —S21 and C2d] B is 
obtained uniquely from c;o; similarly, K corresponds to 
Cdo- 

5. HOT AXISYMMETRIC DISC 

Before we look at the non-axisymmetric case we wish 
to examine the effect of an initially hot stellar disc on 
the measurement of the OC. As described above we use a 
background axisymmetric potential consistent with a flat 
rotation curve. If stars moved on circular orbits and we 
assumed a solar radius of 7.8 kpc and a circular velocity 
of 220 km/s we would measure A = Aq k, 14.1 km/s/kpc 
and B ~ Bq = —Af). This would result in a value for the 
LSR rotation speed fio = A — B = 28.2 km/s/kpc and 
slope of the rotation curve dv^/dr = —{A + B) = 0. 

However, in the case of an initial radial velocity disper- 
sion (T„ = 40 km/s, our simulations produced A « 12.1 
km/s/kpc and B « —14 km/s/kpc. In other words 
B ^ Bq, whereas the hot disk A deviates from its ax- 
isymmetric value by ~ 2 km/s/kpc. This is an inter- 
esting result for even if no perturbation is present we 
already have an error in measuring A. Consequently, 
we have — A — B 26.1 km/s/kpc and dv^/dr = 
~{A + B) ^ 1.9 km/s/kpc. The inferred slope of the 
velocity rotation curve is no longer zero. 

The reason for the observed decrease in A is that our 
hot axisymmetric disk is described not by eqs. [13] but 
by eqs. [121 Radial oscillations cause stars with home 
radii interior and exterior to ro to be observed at the 
LSR. Because the surface density and radial velocity dis- 
persion are exponentially decreasing functions of radius, 
on average, these stars have a lower azimuthal velocity 
than those with guiding radius of rg. This causes a tail 
in the distribution of v weighted toward negative values 
and is quantified by the asymmetric drift, defined as the 
difference betwee n the LSR and the mean ro tation veloc- 
ity: Va = vq — (jBinnev fc Tremaiii3ll987f ). Obviously, 
the asymmetric drift is related to the radial velocity dis- 
persion; it has been determined empirically that in the 
SN V g ~ al/k with A: = 80 ± 5 km/s (jPehnen fc BimievI 
|1998[ ) (fc here should not be confused with the spiral wave 
number). In our notation = —v so we can rewrite eqs. 
[121 as 



2A, 



2B, 



:2^0 



2Bo 



ro 



ro 



dvg 
dr 

dvg 
dr 



(15) 



We numerically estimated the values of the asymmetric 
drift and its slope and found Vg^ro) = 16.9 km/s and 
dva{r)/dr = —1.9 km/s/kpc, respectively. These values. 
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together with the numericaUy derived Aaxi ~ 12.1 and 
-Baxi w — 14 specified above were found to satisfy eqs. 
[T51 The last two terms on the right side of the second of 
these equations happen to almost completely cancel out, 
resulting in _Baxi ~ ^o- 

The deviations from the cold axisymmetric OC (small 
for B and larger for A) estimated from our hot axisym- 
metric disk are in very good agreement wit h the ana- 
lytical results of iKuijken fc Tremaiii^ ()1991l ) and O&D 
(their eq. 12). The asymmetric drift inferred from Va — 
cr'^/k, given a radial velocity dispersion CTm = 40 km/s, is 
20 km/s ±1.3. Even though this is marginally consistent 
with our numerically calculated value Va = 16.9 km/s, 
here consistency is not to be expected since the constant 
k was determined empirically for the SN and thus may 
reflect non-axisymmetric effects or incomplete mixing. 

Eqs. [TO] were derived only to first order in the Taylor 
expansion of the velocities. Higher order terms might be- 
come important for deep samples. In our simulations we 
calculate the OC from a sample of test particles centered 
&t d = 800 pc from the LSR and an annulus width of 200 
pc. Since our expressions for the OC are only valid in 
the limit of d — > 0, we need to evaluate the OC in several 
distance bins and extrapolate to c? = 0. We performed 
such an extrapolation for the hot axisymmetric disk dis- 
cussed above and found that in the range < d < 1.4 kpc 
the variation in A and B was less than 0.5 km/s/kpc. 
However, when spiral structure is included, the variations 
with sample depth could be significant, especially for the 
cold disk case. Sj7]is devoted to this problem. 

The errors in A and B induced by the axisymmetric 
drift in the hot disk need to be kept in mind when consid- 
ering the hot population simulations with included spiral 
structure perturbation, presented in the next section. 

6. THE EFFECT OF SPIRAL STRUCTURE 

We now describe how the presence of spiral perturba- 
tion causes the OC to deviate from their axisymmetric 
values. In all of our simulation runs we used a two-armed 
spiral density wave. To look for variation of the OC with 
the phase of the SN, 0Oi (see Sj4]for definition of 0o) we 
split the arc between two spiral arms into 50 SN disks. 
We estimate the OC in each of these disks as described 
in EH 

It is interesting to see how pattern speeds placing the 
background stars near resonances affect the local veloc- 
ity field. Possible locations o f the LSR have been pro- 
posed to be near the 4:1 ILR (iQuillen fc Minchevll2005f) 
or near the CR (iLepine et al.ll2003HMishurov et al.ll2002t 
iMishurov fc Zeninal 119991 )7" We examine two particular 
cases: 

(i) an initially cold stellar disk and 

(ii) an initially hot disk (cr„ = 40 km/s). 

To find what errors result in measuring the OC when 
spiral structure perturbs the disk, in the figures below 
we plot the difference between A, B estimated from a 
simulation run with an imposed spiral perturbation and 
A, B resulting from an axisymmetric disk. Of course, 
for a cold disk (circular orbits) A = Aq and B = Bq 
but, as we found in the previous section, an axisymmet- 
ric hot disk yields a substantial decrease in A and a small 
increase in B. Therefore, we calculate the spiral struc- 
ture induced errors in A and B differently, depending on 
whether the disk is initially cold or hot: For case (i) we 



plot A A = A - Ao and for case (ii) A A = A - Aaxi. 
Similar expressions hold for AB. 

First we present our results for the variation of the OC 
with phase angle and pattern speed in a more general 
way. In §i j6.3l and l6.4l we concentrate on the 4:1 ILR and 
CR, respectively. Table [T] summarizes the parameters 
used for each run. 

6.1. Cold stellar disk 

We first investigate the effect of spiral structure per- 
turbing a cold stellar disk. In fig. [3] we plot contours 
of AA, AB, C, K, A{A ~ B), and A{A + B), against 
the LSR phase angle, 00 1 and the pattern speed, 17^. We 
define A{A ± B) = AA ± AB. Each horizontal hue cor- 
responds to a simulation with a different spiral angular 
velocity, ilg. Here and AB are the deviations from 
the axisymmetric values Aq and B^. Darker regions cor- 
respond to larger values. Minimum and maximum values 
are indicated above each panel. The zero contour level 
in each panel is indicated by a solid line. We change the 
LSR angle between the convex arm at </)o = 0° and the 
concave one at i^o = 180° (see fig. [2]). Dashed lines show 
the locations of CR and 4:1 LRs. 

There is another interpretation of the y-axis of each 
panel in fig. [3) Since logarithmic spirals are self-similar 
and the velocity rotation curve is flat, changing the pat- 
tern speed at the same radius is equivalent to changing 
the radius and keeping the spiral angular velocity the 
same. Thus, we note that fig. [3] could be perceived as 
the variation of the OC with respect to radius, r/rp. For 
example, assuming the LSR were located at the CR (or 
r^s = r^o) the vertical axis varies between r = 0.6ro and 
r — 1.4ro. Note, however, that this interpretation of the 
y-axis of fig. [3]does not reflect the slope change in the ex- 
ponential surface density distribution which would result 
with the change of radius. 

AA and AB are maximized midway between the spiral 
arms (top left and right panels in fig. [3|). There is some 
structure in A and B as one approaches the 4:1 ILR 
and contours are skewed toward the convex (small (jjo) or 
concave (large (f)o) arm, respectively. 

The Oort's C is directly related to the radial streaming 
caused by the presence of the spiral arms. Its variation 
with 00 and Qs in the cold disk case is shown in the left 
column, second row panel in fig. [3l Structure appears 
to be symmetric across the diagonals with the most no- 
table exception of the clump at {(j)o,Qs) ~ (120°, O.TJIq). 
An interesting result is that C « at the CR for 
all angles. The variation of C with 0o becomes more 
prominent as one moves from the CR to each of the 
4:1 LRs. That is where minimum and maximum val- 
ues are reached: Cmin = —6.3 km/s/kpc at the point 
(0o,r2s) « (20°, 1.35rio) and Cmax = 6.5 km/s/kpc at 
{(j)o,Qs) ~ (25°, O.Tfio). We also note that all structure 
in this panel seems shifted toward small 0o, or the convex 
arm, by about 5°. 

If the LSR were located near the CR then C would 
not provide any information about the spiral structure, 
regardless of the phase angle. On the contrary, A and 
B show marked variation with galactocentric azimuth 
along the CR. However, if in addition we happened to be 
at (po « 45° or <j>o « 135°, then we would measure values 
for A, B, and C (and even K) close to those expected 
for an axisymmetric galaxy. 
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We discuss how an initially hot stellar disk is affected 
by spiral structure in the next section. 

6.2. Hot stellar disk 

In ij5]we showed that a hot axisymmetric disk already 
causes miss- measurements of A and B. We expect this 
result to be reflected in the case of a hot stellar disk with 
an imposed spiral perturbation. To separate only the 
effect of the perturber on A and B we need to correct 
for the asymmetric drift induced error. In the following 
figures we plot = A — A^xi and AA = B — B^xi as 
discussed at the beginning of §6) 

Fig. [4] shows the effect of a two-armed spiral density 
wave on an initially hot axisymmetric background disk. 
The simulation setup is identical to fig. [3]with the excep- 
tion of an initial radial velocity dispersion (t„ = 40 km/s. 
We found in ij5]that this initial setup caused an asym- 
metric drift of Va = 16.9 km/s when non-axisymmetric 
perturbation was absent. 

For the cold disk, maximum AA and A_B values were 
attained in the central regions of their contour plots (top 
panels of fig. [S]). With the increase of velocity dispersion 
both of these have split into two islands, with lines of 
symmetry at ris/^^o ~ 0.9 (fig. [4|). Along this line, and 
even along the CR (see fig. [7]), AA and especially A_B are 
very close to zero. As a direct consequence of the effect 
of the initially hot stellar population, the spiral structure 
induced errors in A and B have decreased by factors of 
^ 3 and ~ 2, respectively. 

Since the contour plots of A A and Ai? in fig. [4] exhibit 
quite similar morphology they cancel out in A(A — B) 
and add up in A{A+B) (bottom two panels in fig. 0]). As 
a result, the A(A — B) contour plot is very similar to its 
cold disk counterpart, whereas the contours of A(A — B) 
resemble those of A A and AS. 

The value of C along the CR line is still close to zero 
and again, the most variation with galactic azimuth takes 
place near the 4:1 LRs. However, similarly to AA and 
A_B, minimum and maximum values have decreased in 
magnitude, by a factor of ~ 2 for C. This was expected 
in view of the fact that with the increase of the stars' 
random motions, the response to spiral structure must 
decrease. The prominent peak seen in the cold C contour 
plot at (0o,f^s) « (120°,0.7r2o) (fig. W) has completely 
disappeared in the hot disk plot (fig. |4|). 

In addition to C w along the CR line found for the 
cold disk, here A A, Ai?, and K are also nearly zero for 
all It is interesting that for the hot disk the symme- 
try line of all contour plots has shifted from the CR to 
ris/rio ~ 0.9. If the Sun is located at, or just inward of, 
the corotation radius then measurements of the OC, us- 
ing a high velocity dispersion population, would provide 
no information about the spiral structure. 

In their observational measurements O&D find that 
C varies approximately linearly with color, B — V, with 
C w km/s for early type and C « —10 km/s for the 
later stars (red giants). They also found similar varia- 
tion with asymmetric drift Va^ where C ~ —10 km/s for 
Va = 18 km/s. Early type stars are younger and have 
small velocity dispersion. On the other hand, high B — V 
suggests older stellar population and thus higher Uu , thus 
larger Va- This contradicts our intuitive expectation, as 
well as our numerical results, that the radial streaming 
induced by spiral structure becomes less important as 



velocity dispersion increases. 

We next look at the variation of the OC with Sun's 
phase angle at a fixed spiral pattern speed. 

6.3. Near the 4-1 ILR with spiral pattern 

Variation of the OC and their combinations as a func- 
tion of LSR phase angle, 0o, at the same spiral pattern 
speed Vis = O.6r2o, or just inside the 4:1 ILR, is shown 
in fig. O Solid and dashed lines represent slices from fig. 
[3] (cold disk) and fig. |4] {au = 40 km/s), respectively. 
The units of the y-axis are km/s/kpc. Deviations from 
axisymmetry (zero y values) decrease with increasing ve- 
locity dispersion in all panels since random motions take 
precedence over the spiral structure perturbation: both 
less structure and decrease of amplitude is observed for 
the hot disk (dashed lines). 

Strong, two-armed spirals have been shown to 
create square-shaped orbits near their 4:1 ILR 
(jContopoulos fc Grosbo]||l986[ ). The orientation of the 
orbits interior to the 4:1 ILR is such that they support 
the spiral structure. On the other hand, at guiding radii 
exterior to the 4:1 ILR orbits of stars fail to support the 
spiral structure. We expect this behavior to be reflected 
in our results for the OC. 

On the outside of the 4:1 ILR, fls — O.lfta, we ob- 
serve a significant change in the variation of the OC for 
the cold disk (solid lines in fig. [6]). However, the OC 
resulting from the hot disk (dashed lines in fig. [6]) are 
almost identical to those in fig. [5l Again, the variations 
present in the cold disk values are reduced by the random 
motions of the stellar population. 

The bottom two panels in figs. [5]and[6]show the third- 
order coefficients in the Fourier expansion of /xj, si^s and 
Q,3 (see These are quite large for the cold disk and 

almost zero for the hot one. As with C, non-zero si^s and 
cj_3 are indicators of non-axisymmetry. 

6.4. At the CR with spiral pattern 

Unlike the prominent variations of the OC with phase 
angle near the 4:1 LRs, at the CR we observe only mild 
effects (fig. [?])• The hot-disk OC values still deviate 
from axisymmetry less than their cold-disk counterparts, 
however the difference is not as pronounced as it is near 
the 4:1 LRs (cf. figs. ElandEl). 

A A and Ai? exhibit similar variation with (j)Q. As a 
result, A(A — B) is very small (bottom left in fig. [7]). 
On the other hand, the inferred slope of the galactic 
velocity rotation curve, dv^/dr — —{A + B), deviates 
from zero except for particular phase angles. It is very 
interesting that C « for all phase angles for both cold 
and hot disks, as mentioned in the discussion of figs. [3] 
and[4l If the Sun is located near the CR, this behavior 
of C would provide no information on the spiral struc- 
ture. The third order Fourier expansion coefficients s;_3 
and c;^3 are very close to zero for both cold and hot disks. 

We ran a simulation with a four-armed spiral structure 
with a pattern speed placing background stars at the 
CR. Our results were similar to the two-armed case. 

In all of our simulations involving a spiral density wave 
perturbation the absolute values of C for the initially 
cold disk were larger than those for the hot disk for 
most values of the phase angle, 0o, and pattern speed 
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0.6 < fis/i^o < 1-4- This shows that it is not possible to 
reproduce the observational measurements of O&D with 
the spiral structure alone. 

7. VARIATION WITH SAMPLE DEPTH 

For all of the simulations with included spiral structure 
discussed above we estimated the OC from a sample at 
an average heliocentric distance d = 0.8 kpc and annulus 
width of 0.2 kpc. We found in SjSlthat for a hot axisym- 
metric disk changing d between 0.2 and 1.4 kpc caused 
almost no variations of A and B. However, when the spi- 
ral density wave perturbation is included this need not be 
the case because there are strong spatial variations in the 
mean velocity field. Following the same approach we es- 
timated the OC in different heliocentric distance bins in 
the range 0.2 < d < 1.4 kpc. The results for fig — O-THq, 
or just outside the 4:1 ILR, and phase angles (j)Q = 68° 
and 00 — 162° are presented in figs. [5]and[51 respectively. 
Solid and dashed lines show the results from simulations 
with initially cold and hot disks, respectively. The y- 
axis units are km/s/kpc. In fig. [5] we observe a large 
variation with sample depth for the cold AA and AB for 
d <~ 500 pc. Due to the similar variation in AA and AB 
mainly the combination A{A + B) is affected, resulting 
in a different value for the inferred slope of the galactic 
rotation curve, depending on the heliocentric distance of 
the sample considered. The cold disk yields much larger 
variations of A and B (about 30% for A) with sample 
depth compared to the hot values (about 10%). 

We also looked at the variations of the OC with average 
distance from the Sun at a phase angle (/iq — 162o, which 
places the SN near to the concave (leading) arm. In fig. 
[9l AA and AB are found to increase with decreasing d, 
which is the opposite to the trend observed in fig. [H The 
result again is a large variation for A{A + B) and smaller 
for A{A — B). However, at this phase angle all other OC 
(resulting from the cold disk), including the third order 
coefhcients, are greatly affected by the change in d. 

We also plotted the variation of the OC with d for a 
simulation with fig = fin and (po — 162° (fig. [TO)) . At 
this pattern speed the usual OC vary with sample depth, 
for both cold and hot disks, by ^ 1 km/s/kpc, however 
the coefficient of the cos 31 Fourier term changes by more 
than ~ 2 km/s/kpc. 

The following argument was found to shed some light 
on the question: Why do the cold OC vary so much with 
the depth of the sample considered whereas the hot ones 
show only mild variations? A problem with the deter- 
mination of the OC addressed by O&D is the effect of a 
non-smooth mean velocity field. Because of local anoma- 
lies in the Galactic potential caused by spiral arms, for 
example, the streaming field would exhibit small scale 
oscillations on top of the underlying smooth field. This 
would result in an oscillating mean velocity field, giving 
rise to significant higher order terms in the Taylor expan- 
sion of the velocities. O&D have simulated this effect in 
their fig. 2 for an axisymmetric streaming velocity. They 
find that for wiggles of wavelength 2 kpc and amplitude 
of only 2% on an otherewise smooth rotation velocity, 
depending on the depth of the stellar sample considered, 
the measured A and B could deviate by as much as 30%. 
A true representation of the smooth background poten- 
tial would be achieved only at distances larger than the 
wavelength of the small scale oscillation. 



To test whether it is a non-smooth streaming field that 
gives rise to the vehement variation of the cold OC with 
sample depth, we plotted the tangential velocity as a 
function of galactic radius in fig. 1111 The x-axis extends 
from (ro — 1.4) to (rp + 1.4) kpc which covers the max- 
imum sample depth considered in figs. [51-ITUl The top, 
middle, and bottom panels represent the spiral struc- 
ture induced wiggles in the mean tangential velocity field, 
v^{r), corresponding to the cold disk values of the OC 
given in figs. [51 [51 and [TUl respectively. We fitted a line 
to each of the plots in fig. [Tl] (dashed lines) with the 
slopes indicated by "s" in the figure. The position of the 
Sun is given by the star symbol at rp = 7.8 kpc. 

As discussed by O&D, estimating the OC when a non- 
smooth streaming field is present would give different 
results for different sample depths. Indeed, considering 
fig. [51 at d w 0.2 kpc the inferred slope of the rotation 
curve is dv^/dr = —{A + B) ^ 3.5, it decreases quickly 
to « —3.5 at d ~ 0.6 kpc and remains approximately 
constant for greater distances. Now looking at the top 
panel of fig. [Tllwe see that the same behavior is apparent 
here. The slope is positive if a sample of stars close 
to the Sun is considered. As 2d becomes larger than 
the wavelength of the velocity curve oscillation (about 
1 kpc), the slope of the smoother background velocity 
field is recovered: s = —3.6 km/s/kpc. For the same 
plot we fitted a line to V(j)[r) of stars with |r — ro| < 
0.8 kpc and found a value of s — —4, consistent with 
the maximum achieved in fig. [8] at d « 0.8 kpc (panel 
showing A{A + B)). 

The middle panel of fig. [TI] corresponds to the cold- 
disk OC values of fig. [H Here the phase angle is 0o = 
162° which causes the spiral structure induced wiggle 
in the velocity curve to have the opposite slope (in the 
range |r — ro| < 1.4) to the top panel of fig. [TI] Again, 
we find a very good agreement with the inferred value 
for the slope given by dv^/dr = —{A + B) in fig. [9] and 
the line fit to shown in fig. [Til At small distances we 
find s < 0; v^{r) is approximately flat for 0.4 < d < 0.6 
kpc; finally, at an average d = 1.4 kpc from the Sun both 
methods estimate s ~ 1.5 km/s/kpc. 

Lastly, the bottom panel of fig. [TT] shows ^^(r) for 
the cold disk at the corotation radius and a phase angle 
00 = 162°. Here v^{r) appears much smoother compared 
to the simulation with fig = O.7r2o- As in the two cases 
above, the variation in ^^(r) is accurately refiected by 
A{A + B) in the corresponding fig. [TOl Note that as in 
the middle panel (same phase angle), the induced slope 
in the region |r — ro| < 1.4 is positive. 

It is important to realize that the slope of the "large" 
scale velocity field we measure for distances greater than 
~ 1 kpc is induced by the asymmetric drift at that par- 
ticular phase angle and is just a larger scale wiggle on 
the otherwise smooth, mean ^^(r). By asymmetric drift 
here we refer to the local, spiral structure induced radial 
streaming, causing stars exterior or interior to ro to be 
seen at rp . Unlike the asymmetric drift of a hot axisym- 
metric disk (see fj5]) , when a non-axisymmetric perturba- 
tion is present Va and its gradient depend not only on the 
galactic radius but also on the azimuth. Since we start 
with initially cold orbits the strong forcing of the spi- 
rals near a resonance (4:1 ILR or CR in this case) could 
cause both positive and negative values for the derivative 
of the asymmetric drift, dva/dr, depending on the phase 
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(jjQ. Consequently, the mean v^{r) could be declining or 
rising as shown in fig. [TlJ In fact, even Va < values 
are possible in a cold galaxy as is the case of the top 
panel of fig. [TT] (since Vcf,{ro) > vq = 220 km/s). This 
was generally found to occur at locations closer to the 
convex arm for both = 0.7f2o and = f2o- In real- 
ity a negative asymmetric drift has never been observed. 
A possible reason for this could be our phase with re- 
spect to the spiral arms. Another explanation could be 
the fact that stars are actually born with some velocity 
dispersion. We simulated a disk with an initial cr„ — 10 
km/s and found that Wq < was still possible. The spi- 
ral structure induced velocity curve wiggles at 0o = 162° 
shown in the middle [Qg = O.T^Iq) and bottom (fls = ^o) 
panels of fig. [Tl] correspond to positive asymmetric drift 
and a phase angle lagging the concave spiral arm by 18°, 
both consistent with what has been proposed to be the 
case for the SN. 

Examining figs. [8l-IT0l we find the largest variation of 
C with sample depth for (po = 162° and ilg = 0.71^0 
(fig. O. We expect this behavior to arise from spiral 
structure induced oscillations in the first two terms on 
the right side of the third of eqs. [10] (assuming the third 
one can be neglected) , similarly to the variations in A and 
B being the result of an oscillating local average velocity 
field, V(j,{r). The interplay of these terms at different 
average distances from the Sun was found to account for 
the variation of C with sample depth at this location. 

Due to the strong variations of the cold-disk OC with 
sample depth, d, we cannot extrapolate to d = as we 
did for the hot axisymmetric disk. Instead, observational 
measurements of the OC should be compared to non- 
axisymmetric models at the same distance from the Sun 
and for the same velocity dispersions. 

We conclude that low velocity dispersion stars (bluer 
samples) would be good tracers for the local variation of 
the mean velocity field. For the hot disk the dependence 
of the OC on sample depth is usually not significant. The 
results of this section show that we cannot neglect the 
variations of the OC with sample depth for cold popula- 
tions. It is remarkable, however, that even when higher 
order terms are present, the cold-disk A and B provide 
a good estimate for the slope of the mean velocity field. 

8. THE OORT CONSTANTS FROM OBSERVATION 

It is interesting to compare our results to the observa- 
tionally deduced values for the OC. From t he th e n avail - 
able proper motions and radial velocities, lOortI ()1927bl ) 
measured A « 19 and B « —24 km/s/kpc. Since then 
many measureme nts have b een done, mainly of A and 
B. The review of iKerr fc LVn dcn-BcU (1986) yields the 
average va lues A = 14.5 ± 1.3 and B = -12.0 ± 2.8 
km/s/kpc. iKuiiken fc Tremaind ()1991[ ) reviewed the ob- 
servational constraints on C and K and concluded that 
C = 0.6 ± 1.1 and K = -0.35 ± 0.5 km/s/kpc. Thus 
these are zero within the errors as in the case of an ax- 
isymmetric Galaxy. 

The release of the Hipparcos catalog (ESA 1997) made 
it possible to re-analyze the local vel ocity field. In their 
work with the Hipparcos Cepheids, iFeast fc Whiteloc3 
(fl99l found A = 14.8 ± 0.8 and B = -12.4 ± 
0.6 km/s/kpc. These va lues are very close to the 
IKerr fc Lvnden-Belll (fl98l averages stated above but 
have smaller errors. iMignardI (|2000f ) derived A = 11.0 ± 



1.0 and B = —13.2 ±0.5 km/s/kpc for early type dwarfs, 
and A = 14.5 ± 1.0 and B = -11.5 ± 1.0 km/s/kpc for 
the distant giants. This increase in \A\ and \B\ with av- 
erage age is consistent with our results for the case of a 
spiral perturbation. At particular phase angles (mostly 
near the spiral arms), figs. EHT] show larger values for 
h otter popula t ions. 

iTorra et al.l (|2000| ) analyzed a large sample of nearby 
O and B stars from the Hipparcos catalog which they 
complemented wit h distance estimates from Stromgren 
photometry ( Hauck fc Mermi lliod 1998.) and radial ve- 
locities. They analyzed the local velocity field in the 
context of the Gould Belt. The latter object is a ring- 
like structure inclined at 16 — 22° with respect to the 
Galactic plane, extending to about 600 pc from the Sun. 
Their inferred values for stars in the range 0.6 < d < 2.0 
kpc are A = 13.0 ± 0.7, B = -12.1 ± 0.7, C = 0.5 ± 0.8, 
and K = —2.9 ± 0.6 km/s / kpc. Using O and B Hippar- 
cos stars, iLindblad et all ([l997h found C = 0.8 ± 1.1 
and K = —1.1 ± 0.8 km/s/kpc for r > 0.6 kpc (or 
outsid e the Gould belt) , in ag reements with lTorra et al.l 
(|2000[ ). IComeron etldj (|1994| ). using B6-A0 stars in the 
range d < 1.5 kpc. found K ss — 1 km/s/kpc in agree- 
ment with lTorra et all (|2000D and lLindblad et all ^1997). 
However, their inferred value for C is clearly negative: 
C — —8.8 ± 1.1 km/s/kpc, and is consistent with the 
"phase mixing" corrected C derived by OfcD. 

We could not account for the last result (large negative 
C from blue stellar samples) by perturbing a stellar disk 
with spiral structure. The rest of the above observational 
measurements of the OC are easily reproduced with the 
influence of a spiral perturbation for a variety of pattern 
speeds and phase angles (cf. figs. EHZ]). 

However, all of the above results could suffer from the 
"mode mixing" problem, pointed out by O&D. They ar- 
gue that non-axisymmetry is the dominant origin of the 
observed differences of the OC between stellar samples. 
If the reason was non-equilibrium effects then one would 
expect a somewhat erratic behavior instead of the clean 
trends seen in OfcD's fig. 9, left-hand panel. The au- 
thors concluded that sample depth is also unlikely to ac- 
count for these changes. We agree with this conclusion 
for the case of the high velocity dispersion population. 
However, as we showed in the previous section, spiral 
structure propagating in a cold disk could cause marked 
variations in all OC as a function of sample depth, d. 

It should be kept in mind that the various refer- 
ences presented here could be subject to various se- 
lection effects. For examp le, the results presented by 
IKerr fc Lvnden-Belll (fl98l are mostly based on North- 
ern hemisphere proper motion and radial velocity sur- 
veys. Th is causes a bias in the OC towards smaller \A\ 
and \B\ ()011ing fc Merrifieldl [l998l ) . In addition, before 
the release of the Hipparcos catalog Oort's B was ob- 
tained in a non-inertial reference frame and should be 
viewed with caution. 

Lastly, the solar motion is another source of error in 
estimating the OC. The motion of the Sun with re- 
spect to the mean velocity vector of a sample of stars 
leads to cos / and sin I terms for the mean radial velocity 
and proper motions as a function of Galactic longitude. 
These n = 1 terms (see §4.ip have been fit to our numeri- 
cally generated stellar sample, so the mean velocity field 
for each simulated measurement has been removed be- 
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fore the simulated OC measurement. One problem con- 
fronting observational surveys is that the mean velocity 
vector with respect to the Sun may in fact depend on 
the stellar sample. By removing the n = 1 terms in our 
simulated measurements we have bypassed this problem, 
however our simulated sample is not affected by distance 
biases. As pointed out by O&D an additional problem is 
that variations in the mean distance of the sample with 
Galactic longitude can lead to aliasing or " mode mixing" 
as well as difficulty in measuring the mean velocity vector 
of the sample. 

9. SUMMARY AND DISCUSSION 

The importance of the Oort constants is in their simple 
relation to the galactic potential, in the case of vanish- 
ing random motions. The circular orbits have velocity 
rd^/dr = = {Aq — BqY . If we could find the causes 
for the deviations, AA and AS, from the "true" values, 
A and B, then we would provide a direct constraint on 
the Galactic potential. 

In this paper we have investigated the effect of spiral 
structure on the measurements of the OC A, B, C, and 
K. We have performed test-particle simulations with 
an initially cold or hot stellar disk and an imposed two- 
armed spiral density wave perturbation. Variations of 
simulated measurements of the OC with pattern speed, 
galactic azimuth, and sample depth have been explored. 

We find that systematic errors due to the presence of 
spiral structure are likely to affect the measurements of 
the Oort constants. Moderate strength spiral structure 
causes errors of order 5 km/s in A and B. Axisymmet- 
ric, high velocity dispersion disks also yield errors in the 
simulated measurements of Oort's A (as a result of the 
asymmetric drift), but they are smaller, of order 2 km/s 
(see fJH). 

If the Sun is located near the corotation resonance 
then, regardless of the phase angle, we would measure 
C w and thus get no information about the spiral struc- 
ture. If, in addition, we constrain (j)Q 45°, or 135°, all 
OC would have values within the current measurement 
error. As the Lindblad resonances are neared, C varies 
strongly with phase angle and has an increased value. 

In |j7]we investigated the effect of sample depth on the 



Amaral, L. H., & Lepine, J. R. D. 1997, MNRAS, 286, 885 
Binney, J., & Tremaine, S. 1987, Galactic Dynamics, Princeton 

University Press, Princeton, NJ 
Comeron, F., Torra, J., & Gomez, A. E. 1994, A&A, 286, 789 
Contopoulos, G., & Grosbol, P. 1986, A&A, 155, 11 
Guddeford, P., & Binney J. 1994, MNRAS, 266, 273 
Dchnen, W., & Binney, J. J. 1998, MNRAS, 298, 387 
Drimmel, R., & Spergel, D.N. 2001, ApJ, 556, 181 
Feast, M., & Whitelock, P. 1997, MNRAS, 291, 683 
Hauck, B., & Mermilliod, M. 1998, A&AS, 129, 431 
Kerr, F. J., & Lynden-Bell, D. 1986, MNRAS, 221, 1023 
Kuijken, K., & Tremaine, S. 1991, Dynamics of Disc Galaxies, 71 
Kuijken, K., & Tremaine, S. 1994, ApJ, 421, 178 
Lepine, J. R. D., Acharova, I. A., & Mishurov, Y. N. 2003, ApJ, 

589, 210 

Lewis, J. R., & Freeman, K. C. 1989, AJ, 97, 139 
Lin, C. C., Yuan, G., Shu, & Frank H. 1969, ApJ, 155, 721 
Lindblad, P. A. B., Kristen, H., Joersaeter, S., & Hoegbom, J. 1997, 
A&A, 317, 36 



measurements of the OC. We found that for an initially 
cold disk, spiral structure can give rise to marked vari- 
ations of all OC with the average heliocentric distance 
of the sample. For A and B this results from an oscil- 
lating average v^{r), whereas C is sensitive to variations 
in the radial velocity and its gradient. We find that for 
phase angles placing the Sun near the convex arm nega- 
tive values for the axisymmetric drift are possible for low 
velocity dispersion stars. Unlike the fT„ induced asym- 
metric drift in a hot axisymmetric disk, in a cold disk it 
is the spiral structure that causes it. It is possible the 
reason for never observing Va < could be our proximity 
to the concave spiral arm. 

Our result for the variation of C with velocity disper- 
sion is distinctly different than the observational mea- 
surement of C done by O&D. Whereas in our simu- 
lations the absolute value of C decreases with increas- 
ing velocity dispersion, the opposite trend is observed in 
fig. 6 by O&D. In the same paper the authors suggest 
that it is possible that the Galactic bar is responsible 
for this behavior. Since orbits change orientation at the 
2:1 OLR of the bar, at this location the bar is expected 
to most strongly distort the local kinematics. As dis- 
cussed by O&D, stars from inside the OLR should pro- 
duce C < 0. If the Sun's location is just outside the OLR 
then we would sample more stars giving rise to negative 
C with increasing sample depth. To test this possibility, 
O&D analyzed only stars brighter than Vr = 10.5 which 
yielded a decrease in the "phase mixing" corrected value 
of C from ^ —10 to ~ —7.5 km/s. Another way to try 
to explain it is to consider the combined effect of spiral 
structure and bar perturbations. Other possibilities are 
the effec ts from minor mergers, or those due to a triax- 
ial halo (jKuiiken fc Tremainelll994( l. Future work should 
aim at understanding this unusual behavior for C. 
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Fig. 1. — A diagram showing the geometry used in the derivation of the OC (see The radial and transverse velocities vectors as 
seen from the Sun are indicated by and d; , respectively. The galactocentric velocity vectors are omitted for clarity. 




Fig. 2. — For a spiral density wave moving with Qs = f2o the corotation (solid) and 4:1 LRs (dashed) circles are shown. The phase angle 
of the SN is given by 0Oi the angle between the galactocentric lines passing through y=0 and the star symbol (representing the Sun). We 
treat 0o as a free parameter in our simulations. In an inertial reference frame the rotation of both spirals and galaxy is in the clockwise 
direction. rf>o increases in the same direction. 
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Fig. 3. — Contours show the spiral structure induced errors in the values of the OC and combinations thereof, as a function of the Sun's 
phase angle, </)o, and the pattern speed, f2s/f7o. A two-armed spiral perturbation with strength = 1.0 is imposed on an initially cold 
stellar disk. The phase angle varies between the convex spiral arm at </<o = 0° , and the concave one at 0o = 180° . Minimum and maximum 
values are given above each panel. Darker regions correspond to larger values. The zero contour level in each panel is indicated by a solid 
line. Dashed lines show the occurrence of resonances. It is interesting to note that C is nearly zero at th e corotation resonance for all 
angles. We note that the y-axis can also be interpreted as a variation of the galactic radius r/ro (see TO. 11 1. Spiral structure parameters 
used in our simulations are found in table [T] 
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Fig. 4. — Same as Fig. [3] but for an initially hot stellar disk. The initial radial velocity dispersion is Cu = 40 km/s. The asymmetric 
drift induced errors in A and B are subtracted so that the plots reflect the effect of the spiral structure only (see beginning of The 
hot population causes the central regions of the cold AA and AB to split into two islands with a line of symmetry along f2s/f2o ~ 0.9. 
In addition to C ~ along the corotation resonance found for the cold disk, here A A, AS, and K are also nearly zero for all 0o- It is 
interesting that for the hot disk the symmetry line of all contour plots has shifted from the CR to Qs/Hq ~ 0.9. If the Sun is located at, or 
just inward of, the corotation radius then measurements of the Oort constants, using a high velocity dispersion population, would provide 
no information about the spiral structure. Deviations from the axisymmetric values of all constants is decreased as a direct consequence of 
the initially hot stellar population. 
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Fig. 5. — Variation of the OC and combinations as a function of the LSR phase angle, <^o, at the same spiral pattern speed, Qs = 0.6f2o 
(just inside the 4:1 ILR). Solid and dashed lines represent slices from fig. |3](cold disk) and fig. |4]((t„ = 40 km/s), respectively. The units 
of the y-axis are km/s/kpc. Deviations from axisymmetry decrease with increasing velocity dispersion in all panels since random motions 
take precedence over the spiral structure perturbation. 
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Fig. 6.— 



Same as fig. [5]except that Qs = 0.7f2o (just outside the 4:1 ILR). 
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Fig. 7. — Same as figs. [5] and [6] but with Qs = f^o (CR). We observe almost the same variation in AA and AB. This gives rise to 
AA — B but large slope of the galactic rotation curve, dv^/dr = —(A + B). It is very interesting that C ~ for all phase angles for 
both cold and hot disks. If LSR is located at or near the corotation resonance l ILepine et al.ll2003l : iMishurov et al.ll200l: IMishurov fc Zeninal 
119991) then C will provide no information about the spiral structure. 
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Fig. 8. — Variation of the OC with sample depth for Qs = 0.7r2o (just outside the 4:1 ILR) and a phase angle 0o = 68° (that is, 
the Sun is 68° ahead of the convex arm). Solid and dashed lines show the results from simulations with an initially cold and hot disks, 
respectively. The x-axis shows the change of the average heliocentric distance of the sample of stars used to estimate the OC. The y-axis 
units are km/s/kpc. Note that in all previous figures the Fourier coefficients presented were estimated at an average distance of d = 0.8 
kpc. At large d the inferred slope of the rotation curve, dv^/dr = —{A + B), is negative whereas for small d it is positive. This is in very 
good agree ment with the shape of the actual wiggle in the average ^^(r) induced by the spiral structure at this location, shown in the first 
panel of fig. 1111 
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Fig. 9. — Same as fig. [8]but for a phase angle </>o = 162°. Here the inferred rotation curve has the opposite behavior to fig. |8] it is rising 
for large d and declining if closer samples are considered. Again, this is in very good agreement with the spiral structure induced wiggles 
in the initially fiat i'i^(r), shown in the second panel of fig. 1111 
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Fig. 10. — Same as figs. |8] and [9] but for a pattern speed Qs = f^o (corotation resonance). The phase angle is <l>o = 162°, same as in fig. 
[9] At this pattern speed the variation of the usual OC with sample depth for both cold and hot disks is only of the order of 1 km/s/kpc, 
however the coefficient of the cosSi Fourier term varies by more than ~ 2 km/s/kpc. 



19 



?35 



If) 



> 



CO 



> 



> 



s = -3.6, 00=68 , U^=Q.1Q,_ 




220 

215 
210 



: s = 4.0, 9!)o=1 6 


7°, n^=nn 







[kpc] 



Fig. 11. — Wiggles in the average tangential velocity v^{r), induced by a two-armed spiral structure propagating in a cold stellar disk. 
The phase angle and pattern speed are indicated in each panel. A linear fit is performed (dashed lines) with the slope values given by "s". 
The Sun's position is at rg = 7 .8 kpc, indicated by the star symbol. The top, middle, and bottom panels correspond to the simulations 
giving rise to figs. |8]|9l and 1101 respectively. 
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TABLE 1 

Parameters describing simulations 





Figures 


4>o 




cTu [km/s] 


d [kpc] 


3 
4 
5 
6 
7 
8 
9 
I 


solid/dashed) 

solid/dashed) 

solid/dashed) 

solid/dashed) 

solid/dashed) 

Jfsolid / dashed) 

ijtop/mid/bottom) 


- 180° 
- 180° 
- 180° 
- 180° 
- 180° 

68° 

162° 

162° 
68°/162°/162° 


0.6 - 1.4 
0.6 - 1.4 

0.6 

0.7 

1.0 

0.7 

0.7 

1.0 

0.7/0.7/1.0 




40 
0/40 
0/40 
0/40 
0/40 
0/40 
0/40 



0.8 
0.8 
0.8 

0.8 

0.8 
0.2 - 1.4 
0.2 - 1.4 
0.2 - 1.4 



Note. — Spiral pattern parameters corresponding to the simulations shown in the Figures. For all simulations there are three parameters 
that we do not change: (i) the spiral perturbation strength is ts = —0.01, given in units of v^, the velocity of a star in a circular orbit at 
ro; (ii) the parameter a = —6, sets the pitch angle of the spiral arms, p, as mcot{p) = a; and (iii) m = 2 since we only consider two-armed 
spirals. 0o is the angle between the galactocentric lines passing through the solar neighborhood and the intersection of the convex spiral 
arm and the circle ro; (j>o = on the line crossing the convex arm. The pattern speed, f!a is in units of = vo/rQ. The radial velocity 
dispersions is given by ctu- Finally, the average distance of the "heliocentric" bins is indicated by d. 



